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Abstract. 

With an emphasis in detecting Earth-like planets set forth by the 2010 Decadal 
Survey and in searching for planets around M dwarfs set forth by the 2008 Exoplanet 
Task Force, radial velocity surveys with infrared echelle spectrometers will have a sig- 
nificant impact on future exoplanet studies. Here, we present the results of an infrared 
radial velocity survey of a sample of 14 late-M dwarfs with the NIRSPEC echelle 
spectrometer on the Keck II telescope. Using telluric lines for wavelength calibration, 
we are able to achieve measurement precisions of 150-300 m/s over a year-long time- 
frame. While we require more RV epochs to determine whether most of our stars have 
planetary-mass companions, we have placed upper limits of 5-10 My on the masses of 
planets around a sub-set of our sample. We have also determined the rotational veloc- 
ities for all the stars in the sample and offer our multi-order, high-resolution spectra 
over 2.0 to 2.4 jjm to the modeling community to better understand the atmospheres of 
late-M dwarfs. 



1. Introduction 

Optical radial velocity surveys have resulted in the detection of hundreds of exoplanets 
with a wide variety of characteristics including many surprises. Now, with measure- 
ment accuracies consistently at the 1 m/s level with improvements in instrumentation 
expected to push that down to 0.1 m/s (Megevand et al. 2010), optical RV surveys 
are capable of detecting super-Earth mass planets around some of our nearest M dwarf 
neighbors. M dwarfs play a critical role in both of these scientific pursuits. The two 
competing planet formation paradigms, core accretion and gravitational instability, pre- 
dict very different planet frequencies for the smallest stars; core accretion models pre- 
dict relatively few planets (1%; Kennedy & Kenyon 2008) while gravitational insta- 
bility models (e.g. Boss et al. 2006) predict frequencies similar to those of solar-type 
stars, perhaps 10% or more. Despite M dwarfs contributing the majority of the planets 
with masses less than 10 M ffi sini, they have been neglected in many of the original large 
scale RV programs since they are so optically faint (V > 14) which degrades the S/N 
of a given RV measurement and thus reduces the measurement accuracy. Finally, M 
dwarfs are, by far, the most abundant spectral type in the stellar neighborhood — 73% 
of all stars within 10 pc are M dwarfs (Henry et al. 2006). 
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One solution to the faintness problem is to move the observations to the near- 
infrared part of the spectrum where M dwarfs are considerably brighter than solar type 
stars. In addition, the contrast between the photosphere of the star and the starspots 
common among M dwarfs is smaller. These co-rotating features are a problematic 
source of noise when looking for the periodic RV signal produced by an orbiting planet. 
The periodicity of activity on the photosphere has been especially problematic when 
searching for planets around young stars (Setiawan et al. 2008; Huelamo et al. 2008; 
Prato et al. 2008). 



2. Observations and Analysis 

We have observed a small sample of 14 late M dwarfs (see Tabled) with the NIRSPEC 
spectrometer on the W. H. Keck II telescope (Mc Lean et al. 1998). In total, the 
observation runs span over three years with some of the targets observed 2006/2007 and 
others observed in 2009. With this instrument we collected spectra in the NIRSPEC-7 
(2.22 /mi) passband in echelle mode with the 3 pixel slit (0.432"), at an echelle angle of 
62.65 and grating angle of 35.50. This configuration results in a wavelength coverage 
of 1.99 to 2.39 lira and a spectral resolution of approximately R=24000. While the 
echelle data has seven full spectral orders, for the radial velocity measurements, we 
focus on the 33rd order covering the spectral range of 2.28-2.32 fim. This wavelength 
region is ideal for this purpose since it is permeated by telluric absorption features as 
well as the CO bandhead prominent in M dwarf spectra. 

Table 1 . Sample of Late M dwarfs 



Target 


RA [2000] 


Dec 


[2000] 


Spectral Type 


K. s 


2M0027+22 


00 


27 


55.9 


+22 


19 


32.8 


M8.0 


9.57 


2M0 140+27 


01 


40 


02.6 


+27 


01 


50.6 


M8.5 


11.43 


2M0253+16 


02 


53 


00.9 


+ 16 


52 


53.3 


M6.5 


7.59 


2M0253+27 


02 


53 


20.3 


+27 


13 


33.2 


M8.0 


11.48 


2M0320+18 


03 


20 


59.7 


+ 18 


54 


23.3 


M9.0 


10.64 


2M1546+37 


15 


46 


05.4 


+37 


49 


45.8 


M7.0 


11.41 


2M1835+32 


18 


35 


37.9 


+32 


59 


54.6 


M8.5 


9.17 


2M1707+64 


17 


07 


18.3 


+64 


39 


33.1 


M9.0 


11.38 


2M1757+70 


17 


57 


15.4 


+70 


42 


01.2 


M7.5 


10.40 


GJ 752 B 


19 


16 


57.6 


+05 


09 


01.6 


M8.0 


8.77 


2M2052-23 


20 


52 


08.6 


-23 


18 


09.6 


M6.5 


11.29 


2M2306-05 


23 


06 


29.3 


-05 


02 


28.6 


M7.5 


10.30 


2M2313+21 


23 


13 


47.3 


+21 


17 


29.4 


M6.0 


10.44 


2M2235+18 


22 


35 


49.1 


+ 18 


40 


29.9 


M7.0 


11.37 


2M2349+12 


23 


49 


49.0 


+ 12 


24 


39.0 


M8.0 


11.56 



The details of the data reduction and spectral extraction are given in Bailey et al. 
(2011 in preparation) and the full description of the results of this survey are given in 
Tanner et al. (201 1 in preparation). Here, we summarize the data analysis and results 
of the survey. After flat fielding and dark correction, the 1-D spectra are extracted 
from the NIRSPEC images using an optimized method first proposed by Home (1986). 
The wavelength calibration and initial values of the properties of the spectrograph are 
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derived from a set of A star calibrators observed over the same night as the science 
targets. A stars are good calibrators in this instance since the spectrum of the star is 
featureless in the wavelength region used to estimate the radial velocities. For some of 
the spectra, we removed bad pixels by hand by replacing the significantly deviant pixels 
with the median value of the surrounding continuum. 



atmospheric transmission spectrum 




n observerd spectrum and best fit model (SNR: 147) 



CO 



E 




residuals (mean: 1.57%) 



2.290 2.295 2.300 2.305 

Wavelength (microns) 

Figure 1. Plot of one epoch of the telluric (top), synthetic model (middle) and 
combined, modeled (bottom-red) and observed (bottom-black) spectra used to deter- 
mine radial velocities for 2M 1757 +70. The bottom, dotted line shows the residuals 
after the subtraction of the modeled and observed spectra which have an rms of 
-1.6%. 

The observed spectrum is the combination of the intrinsic spectrum of the star and 
the telluric spectrum convolved with the instrument point spread function (PSF). We 
utilize a telluric spectrum derived from a solar spectrum (Livingston & Wallace 1991) 
and NextGen model stellar atmospheres of the late M stars derived as a function of 
effective temperature and surface gravity (Hauschildt et al. 1999). For the A stars, 
we assume a flat, featureless intrinsic stellar spectrum. Figure Q] shows the best fit to 
the data for the June 27, 2007 epoch. Also, plotted are a model atmospheric transmis- 
sion spectrum and the theoretical spectrum of an M8 dwarf. The bet fit to the data is 
determined through x 1 minimization. 



3. Results 

Out of the 14 late-M dwarfs in our science sample, five of them have multiple obser- 
vations allowing us to determine the rms of their nightly radial velocity measurements. 
For instance, the late-type M dwarf, GJ 752 B has six epochs of RV measurements over 
the span of five months (see Figure [2]). These radial velocity measurements have an 
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rms of 250 m/s. This is worse than the 50 m/s we have been getting for our early-M 
dwarf calibration stars (Bailey et al. 2010, in preparation) and the 100 m/s rms that 
we have estimated assuming a perfect instrument and no intrinsic stellar jitter (But- 
ler et al. 1996). This could be due to the lower S/N of these observations or needed 
improvements in the synthetic model for the atmosphere of the M dwarf. 
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Figure 2. Plot of the radial velocity measurements collected for GJ 752 B over a 
five months. These observations have an RMS of 250 m/s. 



3.1. Mass Limits Determined from RV Data 

Two of the stars, GJ 752B and 2M 1757+70, have six and seven epochs of RV data, 
respectively, over a time scale of more than a year. This allows us to constrain on the 
mass of any undetected planets orbiting these stars through Monte Carlo simulations. 
Since this survey is most sensitive to short period companions, we set limits only at 
representative orbital periods, namely 10 and 100 days. For each star, a set of 10,000 
orbits are generated with random inclinations, phases, and companion masses ranging 
from ~ 1 to 30 Jupiter masses. Conservative detection limits are then set by determining 
the companion mass that induces a RV dispersion that is 2<r greater than the expected 
dispersion for the star's vsim 99% of the time. With these assumptions, we can rule out 
planets around these two stars with masses of more than ~6 and 10 My at periods of 10 
and 100 days, respectively. 

3.2. Rotational Velocities of Late-M dwarfs 

Because the rotational velocity is one of the free parameters used when modeling the 
observed spectrum, we are able to determine this parameter for all of the stars in our 
sample. Figure [3] plots the rotational velocities estimated from our infrared spectra as 
a function of spectral type. In some cases these values are similar to those estimate 
from optical spectra (Reiners & Basri 2010) with good agreement between the infrared 
and optical rotational velocities for the slow rotators (<20 m/s). However, the two 
fastest rotators in the sample have disparaging infrared and optical rotational velocities. 
Reiners et al. 2009 and other have used starspot models to show that the radial velocity 
jitter of a star is directly related its rotational velocity. Therefore, determining the 
rotational velocities for a large portion of the thousands of M dwarfs within 25 pc and 
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proving a rotational velocity versus radial velocity precision correlation observationally 
would greatly benefit future planet search efforts. 




7 8 
M Spectral Type 



Figure 3. Plot of our estimated rotational velocities as a function of stellar spectral 
type. These results resemble the larger sample of Reiners & Basri (2009) which show 
no correlation. We do have evidence that the precision of the RV measurements 
increases with rotational velocity. 



4. Summary and Recommendations 

With radial velocity precisions of 150-200 m/s for stars of masses of 0.08-0.15 M , in 
general, these observations are sensitive to planets with masses of 5-10 Mj and 10-15 
Mj at periods of 10 and 100 days. Considering that there are a few thousand known 
M dwarfs within 25 parsec and these stars represent a population of stars largely un- 
explored with most planet detection techniques, these observations open the door to 
potentially fruitful new discoveries. While the percentage of early-M dwarfs with gas 
giant planets is currently estimated at 3% (Johnson et al. 2010), we have no information 
on the same value for late-type M dwarfs. With the masses of these stars approaching 
that of the most massive brown dwarfs and reaching a mass ratio of Y compared to 
a Jupiter mass planet, the planet population around these stars could have significant 
implications on the method of planet formation that is dominant in this stellar mass 
regime. 

These spectra can be used for more than just radial velocity measurements. Our 
data represent some of the only high-spectral resolution near-infrared spectra of late-M 
dwarfs published to date. While the radial velocity program uses just the 2.28-2.31 
/im order, our NIRSPEC data cover 1.98-2.39 yum (see Figure |4|. With the discovery 
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Figure 4. Plot of all the observed spectral orders with our NIRSPEC setup of the 
M3.5 dwarf GJ 628. These spectra include both the CO bandhead and Br y. 



of brown dwarfs a few decades ago, modeling the spectral features and, therefore, the 
complex atmospheres of late-M dwarfs became overlooked compared to what's be ac- 
complished for solar-type stars. Our M dwarfs spectra are available to any who wish to 
utilize them for spectral synthesis. None of our spectral show signs of Br y emission 
(2.166 jum) which could be used to aid in age determinations. 

With so few targets in this pilot study, its difficult to determine whether there is any 
observational correlation between radial velocity precision and the rotational velocity 
or spectral type of the star. However, we have proven that we are capable of detecting 
Jupiter-mass planets around these stars which are too faint to be observed efficiently 
with optical echelle instruments. Planned upgrades to the NIRSPEC instrument which 
include adding a methane gas absorption cell to improve wavelength calibration could 
lead to radial velocity precisions of -20-30 m/s (Bean et al. 2009). Such observations 
would be sensitive to Jupiter mass planets in orbits of less than a year around late-M 
dwarfs. 
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